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The use of lattice QCD to understand hadron structure is described, with particular 
emphasis on exploring the role of glue. 



1 Introduction 

o 

This conference on excited states of the nucleon emphasizes the wealth of infor- 
mation about hadron structure accessible empirically through clever experiments 
' and frontier electromagnetic probes. For a theorist, it also highlights the concomi- 

\& . tant need for understanding the remarkably rich structure of hadrons from first 

principles. 

Lattice QCD is the only known way to solve, rather than model, QCD. Given 
' the development of the tools of lattice field theory and advances in computer tech- 

nology that now make Terascale lattice computation affordable, the time is right for 
a major initiative in hadronic physics. Thus, it is particularly noteworthy for partic- 
ipants at this conference at Jefferson Lab that Jefferson Lab, MIT, and the Lattice 
Hadron Physics Collaboration, representing twenty-two theorists at thirteen insti- 
tutions, have undertaken such an initiative and proposed dedicated Lattice QCD 
clusters with a total capacity of 480 Gflops to study hadron structure. The un- 
derstanding and support of the experimental community for such an initiative is 
crucial, and I hope to provide in this talk glimpses of some of the physics that will 
be possible. 

' One of the most striking features of the quark and gluon structure of the nu- 

cleon is the extremely important role of gluons as essential dynamical degrees of 
freedom. In most familiar many-body problems, the bosons whose exchange medi- 
ates the interaction between fermions may be subsumed into an effective two-body 
interaction and hence do not appear directly as essential degrees of freedom. In 
atomic physics, photon exchange is subsumed into a static Coulomb potential and 
one solves the Schrodinger equation in the presence of Coulomb forces. Similarly 
in nuclear physics, to a good approximation nuclear structure may be calculated 
as a many-body problem in the presence of a strong, state-dependent but static 
nucleon-nucleon potential. In both cases, the mass of an atom or nucleus comes 
overwhelmingly from the masses of the fermionic constituents. 

Nucleons are strikingly different. We know experimentally that approximately 
half of the momentum and angular momentum of a nucleon is carried by its gluons. 
Furthermore, the mass of a nucleon would hardly be altered if the masses of its 
quarks decreased from their physical scale of the order of a few MeV to zero, so 
the mass arises almost completely from the gluonic interactions between fermions. 
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Hence, although we normally only think of the fermionic component of the atomic 
or nuclear wave function, in the case of the nucleon it is essential to understand 
the gluon component of the wave function. As I will describe below, lattice QCD 
provides an excellent framework to explore this physics. 

In lattice QCD, an observable is evaluated by defining quark and gluon variables 
on the sites and links of a space-time lattice, writing a Euclidean path integral of 
the generic form 



= Z- 1 f T>(U) t P** tM M-s(u) M - 1 (U)M- 1 (U) 



and evaluating the final integral over gluon link variables U using the Monte Carlo 
method. The link variable is U = e maA ^ x \ the Wilson gluon action is S(U) — 
tt$^p(1 — -^RcTrC/n) where £7 P denotes the product of link variables around 
a single plaquette, and M(U) denotes the discrete Wilson approximation to the 
inverse propagator such that in the continuum limit, M(U) — > m + (j) + ig^L. For 
those unfamiliar with lattice QCD, an elementary introduction at the level required 
for this work is presented in refEJ For our present purposes, it is useful to note that 
evolution in Euclidean time provides a convenient filter to extract the ground state 
from an arbitrary state of specified quantum numbers, since 

e -PH^ = ^ e-P^Cnlpn -> const X 00 + 0( e -/3(®i-Eo)). 



2 Lattice Wave Functions 



The wave function plays a central role in our understanding of many-body systems. 
Familiar examples in nonrelativistic quantum mechanics include the shell model for 
nuclear and atomic systems, the BCS wave function for superconductivity and the 
phases of liquid He, and the Laughlin wave function for the quantum Hall effect. It 
is thus natural to seek analogous understanding of hadronic structure in terms of 
quark and gluon wave functionals, which become wave functions when defined on a 
discrete lattice. 

In nonrelativistic quantum mechanics, we define an N-particle wave function as 
the overlap between a state in which the positions of N particles are specified and 
the state \ip), (xi, . . . , xjvIV') = VK^ii • ■ ■ i x n)- When it is impractical to deal with 
the exact wave function, one may alternatively study a variational wave function 
specified by some set of variational parameters {a±, . . . , a n } and evaluate the inner 
product Atrial (c*i; • • • j ^n)!^)' The "best" approximation is obtained by varying the 
a's to maximize the overlap. A familiar nuclear physics example is approximating a 
light nucleus by a harmonic oscillator wave function, and varying the size parameter 
to obtain the maximum overlap. 

For pure Yang-Mills gauge theory on a lattice, we specify the gauge field A fl (x n ) 
for each space-time direction \x and site x„ , so the lattice wave function is the over- 
lap (.Ao(xi), . . . , A^xn)^) corresponding to the wave functional ^[^(x)] in the 
continuum limit. In QCD, it is convenient to use occupation number representa- 
tion for the quarks, and indicate the lattice sites that contain quarks or antiquarks, 
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Figure 1: Lattice calculation of the meson wave function, ip a , corresponding to the overlap between 
the meson ground state on the left with a configuration containing a q-q pair connected by a string 
of glue on the right. 



so that the wave function is written: (q(xi), . . . , q{x n ),q(y\), . . . , q(y m ), Ao(xi), . . . , 
A4 (x jv) I ip). Whereas it will be straightforward to examine components of the wave 
function involving a small number of quarks and antiquarks at specified positions, 
it will not be practical to keep track of all the components of the gauge field A fl (x n ) 
(which exceed 10 s on a small 16 4 lattice). Hence, in general, we will need to con- 
sider some simple or physical Ansatz for A^(a; n ), and I will give illustrative examples 
below. 

Given the form of a trial wave function 

\q{x x ), . . .,q{x n ),q(yi), . . . ,q(y m ),A (x 1 ), . . .,A 4 (x N )) 

the overlap with the ground state hadron wave functions is calculated on the lattice 
as sketched schematically in Fig. [lj A source with the quantum numbers of the 
desired meson is placed at the leftmost time ti , and propagation to time t near the 
middle of the lattice filters out the ground state. On time slice t, one constructs an 
operator corresponding to the trial wave function. In the case considered in Fig. ^, 
a quark and antiquark are separated by distance y, the cm variable x is integrated 
to project onto zero momentum, and the gluon wave function is the ground state 
vacuum configuration of gluons multiplied by a string of flux created by a Polyakov 
line as shown. Summing over an ensemble of gluon configurations then yields the 
quantum mechanical result for the overlap. 

To elucidate the role of the gluon component of the hadron wave functional, 
it is instructive to display the results for several alternative approximations to the 
gluon field. To highlight the basic ideas, I will show some pedagogical results cal- 
culated previously in the meson sector by Kien Boon TecoH. With the current 
interest in nucleon structure and the resources presently available, analogous calcu- 
lations should certainly be carried out for nucleons. In particular, it is instructive 
to consider the following three definitions of variational wave functions. 
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2.1 Axial gauge or string wave function i/j s 

The conventional gauge-invariant wave function, which corresponds to the Bcthc- 
Salpeter amplitude in the axial gauge A y — 0, is given by 

My) = <o|5(o)re/o ' dzAAz) q ( y )\h) (i) 

where T depends on the Dirac structure of the ground state hadron \h). The ex- 
ponential corresponds to a product of U links as shown in Fig. [I] and the implied 
gluonic component of tp s is that of an infinitely thin string of glue. The square of 
the resulting wave function specifies the probability that a qq pair is separated by 
a distance y and that the glue forms a narrow string connecting them. 

2.2 Coulomb gauge-fixed amplitude %j) c 

Gauge fixing to Coulomb gauge surrounds each quark with the gluons corresponding 
to the static Coulomb field. In the Abelian case we can write down the following 
explicit expression: 

MV) = (0\m^J d3zS ^ g) - X{g) q (y)\h) (2) 

where -Ecstatic is the static Coulomb field of an e + e _ pair separated by distance y. 
Since 

J a- s^tic = J a- V(f> = - y (v • A)<t> 

it is evident that the exponent vanishes in Coulomb gauge. 




Figure 2: Lattice calculation of the meson wave function, ip c , corresponding to the overlap between 
the meson ground state (left) with a vacuum gluon configuration in Coulomb gauge (right). 



For QCD, the implied gluonic component of ip c is obtained by the lattice mea- 
surement shown in Fig. ^ with the gluon field at time slice t fixed to Coulomb 
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Gauge. As discussed later, the Coulomb gauge result has a well-behaved continuum 
limit. 



2.3 Adiabatic Wave Function ipa 

A more physical definition of a gluon wave function is the adiabatic wave function 
obtained by letting the gluonic component be the ground state distribution of gluons 
in the presence of a qq pair at separation y: 



ip a (y) = lim 



(0\q(0)rS y (t,n)g(y)\h) 
W(y,2n)i 



(3) 



where S y (t,n) = u^'*^°' t+n) • u^t+n^t+n) . ^(y^t+n) and w ^ 2n ) is a 
Wilson loop of size y x 2n. When the temporal links are extended (by increasing n) 
far enough, the gluons in the vacuum adjust themselves to the presence of the qq 
pair, forming the ground state of a flux tube connecting two fixed color charges. 




Figure 3: Lattice calculation of the adiabatic meson wave function, ip a , corresponding to the 
overlap between the meson ground state on the left with a configuration containing a static qq 
pair connected by a ground state flux tube on the right. 



3 Density-Density Correlation Functions 

Since the density-density correlation function is a gauge-invariant quantity that 
describes the distribution of quarks in a hadron and reduces to the square of the 
wave function in nonrelativistic quantum mechanics, we have for comparison also 
calculated 

(p(x)p(0)) = {h\u l0 u{x)d lo d{0)\h) (4) 

and its projection {pp) q q onto the qq subspace. Computationally the projection is 
performed by using hard-wall boundary conditions to exclude those paths in which 
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Figure 4: Two representative time histories contributing to the meson density-density correlation 
function. Diagram (a) represents a contribution from the qq sector whereas (b) corresponds to the 
~qqqqqq sector. 



the propagators can cross the time slice t at which (pp) is evaluated. Figure ^ 
shows representative time histories corresponding to the density-density correlation 
function in the qq sector, diagram (a), where only one quark and antiquark propa- 
gator cross time slice t, and a general contribution, diagram (b), containing multiple 
quark-antiquark excitations at time t. 

4 Lattice Calculation of Meson Wave Functions 
and Correlation Functions 

Figures ||-0show the results of calculating the three wave functions and .two density- 
density correlation functions described above in quenched lattice QCD0U. 




Figure 5:rPion ip 3 and ip c for (3 = 5.7 (Teo and Negelea) and f} = 6.0 (extracted from Hecht and 
DeGrando) versus qq separation y. The lines are to guide the eye. 
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Figure 6: Comparison of and density-density correlation functions for the pion as a function 
of the qq separation in units of 0.2 fm. Lines connecting Monte Carlo results are to guide the eye. 

All propagators were calculated with twenty quenched configurations at /3 = 5.7 
on a 16 3 x 16 lattice using the Wilson action for both gluons and quarks and quark 
masses of 40, 95, and 170 MeV. We computed ip s using a string of links and explicitly 
fixing to Coulomb gauge for ip c . T = 75 and T = 7^ are used for tt and p mesons 
respectively. 

Figure [| compares ip s and ip c calculated on the lattice for two different values 
of (3 with a quark mass 170 MeV. The f/'c's agree to within error bars, indicating a 
well-behaved continuum limit. In contrast, ip s decreases with increasing f3 because 
of lattice artifacts that diverge as the lattice spacing approaches zero. As discussed 
below, the small overlap between a thin string of glue and the gluons in the ground 
state hadron accounts for much of the falloff of ij) s with qq separation. 

Figures |^ and show the results for the squares of wave functions and density- 
density correlation functions for tt and p mesons, respectively. All quantities are 
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Figure 7: Comparison of \ip\ 2 and density-density correlation functions for the rho meson as a 
function of the qq separation in units of 0.2 fm. Lines connecting Monte Carlo results are to guide 
the eye. 



normalized to 1 at zero qq separation. We display results for a quark mass of 40 
MeV, where errors are in better control, instead of the extrapolated chiral limit 
which differs inconsequentially. 

The dramatic difference between i/j^t V'c) an d "0s shows that the gluon flux tube 
of the adiabatic wave function produces a much larger overlap with the gluons in the 
ground state hadron than does the Coulomb or string wave function. Beyond 1 fm, 
the adiabatic flux tube is favored by more than an order of magnitude. Physically, 
there is a substantial probability to find a quark-antiquark pair separated by 1 fm. 
connected by a physical adiabatic flux tube, but a very small probability to find 
them joined by an unphysically narrow string of flux. Also note, as mentioned in 
connection with Fig |^, that the string of flux becomes increasingly unfavorable as 
one approaches the continuum limit. 
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The role of multiple quark-antiquark excitations in the density-density correla- 
tion function is evident from the fact that (pp) falls off more slowly than {pp) q q at 
large distances. For both the pion and rho, i\>\ is close to (pp) q q, suggesting that the 
adiabatic flux tube is a reasonable approximation to the full gluon wave function. 

5 Overlap with Variational Wave Functions 

A complementary means of studying the overlap of a trial gluon wave function with 
a hadronic ground state is through the calculation of two-point correlation functions. 
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Figure 8: Ratio of excited state content to ground state content of Atrial as a function of the size 
of the trial function in units of 0.1 fm. 

Let ,P be a source with the desired quantum numbers and denote Atrial) = J^|0). 
For example, could be built from point sources, Gaussian distributed sources in 
Coulomb gauge, or gauge-invariant Wuppertal-smeared sources^ of the form: 

/ \ N 
i>{x) = ^{l + u^[U l tx)8 x ^ x+ , + u}{x-i)5 x ,, x -.]} q(x') 

x' i 

where q(x) is the quark creation operator. The two-point function may be written 

(J(T)jt(O)) = C^K^trial | n)| V £ » T . 

n 

At T = 0, the value of the correlation function is A — cJ2 n KV'triail^)! 2 - The large 
T behavior of the correlation function is c|(V'triai|0)| 2 e _£ ' oT , and extrapolating this 
behavior back to T = yields B — c|(?/;triai|0}| 2 . Hence one may study the overlap of 
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Atrial with the physical ground state by calculating B/A — (Atrial |0)| 2 and optimize 
it with respect to the variational parameters in Atrial 

Figure || shows the ratio of overlap with excited states to the overlap with 
the ground state, (A-B)/B = J2 n >o KV^riail^lVKV'triadO)! 2 for a nucleon cre- 
ated with Wuppertal smeared sources of different rms radii] . The striking fact here 
is that, whereas for point sources the excited state content is 10 4 larger than the 
ground state content, for the optimal smearing the overlap with the ground state, 
B/A = | (Atrial 1 0) | 2 , is approximately 65%. In other words, the sum of the sets 
of gauge-field links generated by the smearing provides a reasonable physical ap- 
proximation to the behavior of the gauge fields in the nucleon. This suggests that 
a systematic study of variational wave functions has the potential to yield even 
larger overlaps and to provide a valuable tool for exploring the nucleon's gluon 
wave function. 



6 Conclusions 



Lattice QCD provides a powerful tool for exploring hadron structure and, in partic- 
ular, the gluon wave functional. The exploratory results for meson wave functions 
shown here display large effects of the gluonic components in three definitions of 
hadron wave functions, and the adiabatic flux tube wave function yielded the most 
physical results. Quantitative measurement of the overlap of a trial wave func- 
tion with the ground state may be obtained from two-point correlation functions, 
providing the opportunity for quantitative study of variational wave functions. 
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